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It has recently been shown that traversable wormholes may be supported by phantom energy. In 
this work phantom wormhole geometries are modelled by matching an interior traversable wormhole 
solution, governed by the equation of state p — ojp with uj < — 1, to an exterior vacuum spacetime 
at a finite junction interface. The stability analysis of these phantom wormholes to linearized 
spherically symmetric perturbations about static equilibrium solutions is carried out. A master 
equation dictating the stability regions is deduced, and by separating the cases of a positive and 
a negative surface energy density, it is found that the respective stable equilibrium configurations 
may be increased by strategically varying the wormhole throat radius. The first model considered, 
in the absence of a thin shell, is that of an asymptotically flat phantom wormhole spacetime. The 
second model constructed is that of an isotropic pressure phantom wormhole, which is of particular 
interest, as the notion of phantom energy is that of a spatially homogeneous cosmic fluid, although 
it may be extended to inhomogeneous spherically symmetric spacetimes. 
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I. INTRODUCTION 

It is now generally accepted that the Universe is under- 
going an accelerated phase of expansion 0,0,1113, where 
the scale factor obeys a > 0. This cosmic acceleration is 
one of the most challenging current problems in cosmol- 
ogy. Several candidates, responsible for this expansion, 
have been proposed in the literature, namely, dark en- 
ergy models, generalizations of the Chaplygin gas, mod- 
ified gravity and scalar-tensor theories, tachyon scalar 
fields and braneworld models, amongst others. The dark 
energy models are parametrized by an equation of state 
given by oj = p/p, where p is the spatially homogeneous 
pressure and p is the dark energy density. For the cosmic 
expansion, a value of oj < —1/3 is required, as dictated 
by the Friedman equation a/a = —4iv(p+p/3). A specific 
exotic form of dark energy, denoted phantom energy, has 
also been proposed, possessing the peculiar property of 
oj < — 1. This parameter range is not excluded by ob- 
servation, and possesses peculiar properties, such as the 
violation of the null energy condition and an infinitely in- 
creasing energy density, resulting in a Big Rip, at which 
point the Universe blows up in a finite time How- 
ever, recent fits to supernovae, CMB and weak gravita- 
tional lensing data indicate that an evolving equation of 
state oj crossing the phantom divide —1, is mildly favored, 
and several models have been proposed in the literature 
HHBSHIUHEEIHEI- In particular, models 
considering a rcdshift dependent equation of state, oj(z), 
provide significantly ameliorated fits to the most recent 
and reliable SN la supernovae Gold dataset |l7j . 

As the phantom energy equation of state, p — tup with 
oj < — 1, violates the null energy condition, p + p < 0, 
the fundamental ingredient to sustain traversable worm- 
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hole [l8L H^. l2(i| , one now has at hand a possible source 
for these exotic spacetimes. In fact, this possibility has 
recently been explored 0, , and it was shown that 
traversable wormholes can be theoretically supported by 
phantom energy. However, a subtlety needs to be pointed 
out, as emphasized in Refs. [2lll22| . The notion of phan- 
tom energy is that of a homogeneously distributed fluid. 
When extended to inhomogeneous spherically symmet- 
ric spacetimes, the pressure appearing in the equation of 
state is now a radial pressure, and the transverse pres- 
sure is then determined via the field equations. In this 
context, it is interesting to note that wormhole solutions 
with an isotropic pressure were found in 22], although 
these geometries are not asymptotically flat. Sushkov, 
in Ref. plj . found wormhole geometries by considering 
specific choices for the distribution of the energy density, 
and in |22|. a complementary approach was traced out, 
by imposing appropriate choices for the form function 
and/or the redshift function, and the stress-energy ten- 
sor components were consequently determined. In Ref. 
|22| it was also shown, using the "volume integral quan- 
tifier" [23l I24I ] , that these geometries can be theoretically 
constructed with infinitesimal amounts of averaged null 
energy condition violating phantom energy. 

It is also of a fundamental importance to investigate 
the stability of these phantom wormhole geometries (It 
is also interesting to note that a stability analysis of a 
specific class of traversable wormholes was carried out 
in Ref. [25j. in a rather different context). As in Ref. 
[22j , we shall model these spacetimes by matching an in- 
terior traversable wormhole geometry with an exterior 
Schwarzschild vacuum solution at a junction interface 
|26t l27l 123 . 129^ . In this work, we analyze the stability 
of these phantom wormholes to linearized perturbations 
around static solutions. Work along these lines was done 
by considering thin-shell Schwarzschild wormholes, using 
the cut-and-paste techniq ue 1301 . It was later shown that 
the inclusion of a charge [31j and of a cosmological con- 
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stant [32| significantly increases the stable equilibrium 
configurations found in Ref . 30] . The advantage of this 
analysis resides in using a parametrization of the stabil- 
ity of equilibrium, so that there is no need to specify a 
surface equation of state. Note that the stability anal- 
ysis of these thin-shell wormholcs to linearized spheri- 
cally symmetric perturbations about static equilibrium 
solutions was carried out by assumin g th at the shells re- 
main transparent under perturbation 33] . This amounts 
to considering specific spacetimes that do not contribute 
with the momentum flux term in the conservation iden- 
tity, which provides the conservation law for the surface 
stress-energy tensor. The inclusion of this term, corre- 
sponding to the discontinuity of the momentum imping- 
ing on the shell, severely complicates the analysis. How- 
ever, we shall follow the approach of Ishak and Lake [33| . 
with the respective inclusion of the momentum flux term, 
and deduce a master equation responsible for dictating 
the stability equilibrium configurations for the specific 
phantom wormhole geometries found in Ref. |22|. We 
shall separate the cases of a positive and a negative sur- 
face energy density, and find that the stability may be 
significantly increased by varying the wormhole throat. 

This paper is outlined in the following manner. In 
Section II, we present solutions of a phantom energy 
traversable wormhole. In Section III, we outline a general 
linearized stability analysis procedure, and deduce a mas- 
ter equation dictating stable equilibrium configurations. 
We then apply this analysis to phantom wormhole ge- 
ometries and determine their respective stability regions. 
Finally in Section IV, we conclude. 



II. PHANTOM ENERGY TRAVERSABLE 
WORMHOLES 

A. Field equations 

The interior wormhole spacetime is given by the fol- 
lowing metric 0] 



ds 2 = _ e 2*(r) df 2 + 



dr 2 



1 — b(r)/r 



(d9 2 + sm 2 6d(j> 2 ), 
(1) 

where and b(r) are arbitrary functions of the radial 
coordinate, r, denoted as the redshift function and the 
form function, respectively |18j . The wormhole throat 
is located at b(ro) = r = tq. For the wormhole to be 
traversable, one must demand that there are no horizons 
present, which are identified as the surfaces with e 2 * — > 0, 
so that $(r) must be finite everywhere. The condition 
1 — fe/r>0is also imposed. The stress-energy tensor 
components are given by (with c = G — 1) 
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where p(r) is the energy density; p r (r) the radial pres- 
sure; andpt(r) the transverse pressure. The conservation 
of the stress-energy tensor, T^. u — 0, provides us with 
the following relationship 



Pr 
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(5) 



A fundamental ingredient of traversable wormholcs 
and phantom energy is the violation of the null energy 
condition (NEC), which is defined as T^k^k" > 0, where 
A: M is any null vector and the the stress-energy ten- 
sor. Note that for phantom energy, governed by the equa- 
tion of state lu = p/p with u> < — 1, one readily verifies 
that the NEC is violated, i.e., p + p < 0. For worm- 
hole spacetimes, consider an orthonormal reference frame 
with k* 1 — (1, ±1, 0, 0), so that we have 
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8^ 



b'r-b 



2 1 



b\ ^ 
r ) r 



(6) 



Thus, using the flaring out condition of the throat, (b — 
b'r)/2b 2 > [la, Hil' an d considering the finite character 
of $(r), we verify that evaluated at the throat the NEC 
is violated, i.e., T^k^k 1 ' < 0. Matter that violates the 
NEC is denoted as exotic matter. 

Note that the notion of phantom energy is that of 
a homogeneously distributed cosmic fluid. However, as 
emphasized in 0, [2^ , it may be extended to inhomo- 
geneous spherically symmetric spacetimes by regarding 
that the pressure in the equation of state p = top is now 
a radial pressure p r . The transverse pressure pt may then 
be determined from the field equation, in particular, from 
Eq. Thus, to find phantom energy traversable worm- 
hole spacetimes, we use the equation of state p r = cup 
with to < — 1, representing phantom energy, and thus 
deduce the following relationship 



*'(r) 



2r 2 (1 - b/r) ' 



(7) 



by taking into account Eq. @ and Eq. ©. 

To model a traversable wormhole, one now considers 
appropriate choices for b(r) and/or <E>(r). Note that this 
is necessary as we only have four equations, namely, Eqs. 
101 -®, and Eq. (0, with five unknown functions of r, 
i.e., p(r), p r (r), pt{r), b(r) and $(r). We shall only con- 
sider form functions of the type b'(r) > 0, as in cosmology 
the phantom energy density is considered positive. Now, 
using the flaring out condition evaluated at the throat 
|l8ll20| . we also have the condition b'(ro) < 1- 

One may construct asymptotically flat spacetimes, in 
which b{r)/r — > and $ — > as r — > oo. However, one 
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may also consider solutions with a cut-off of the stress- 
energy, by matching the interior solution to an exterior 
vacuum spacetime, at a junction interface, a. For sim- 
plicity, in this paper, we shall consider that the exte- 
rior spacetime is the Schwarzschild solution, so that the 
matching occurs at a junction interface, r = a, situated 
outside the event horizon, i.e., a > r& = 2M, in order to 
avoid a black hole solution. 



B. Specific phantom wormhole models 

The physical properties and characteristics of specific 
phantom energy traversable wormhole models were an- 
alyzed in Ref. by considering asymptotically flat 
spacetimes and by imposing an isotropic pressure. Using 
the "volume integral quantifier" it was found that it is 
theoretically possible to construct these geometries with 
vanishing amounts of averaged null energy condition vi- 
olating phantom energy. Specific wormhole dimensions 
and the traversal velocity and time were also deduced 
from the traversability conditions for a particular worm- 
hole geometry. We shall briefly summarize two specific 
phantom wormhole models, found in Ref. [2^, and for 
which we shall further analyze the respective stable equi- 
librium configurations. 



Asymptotically flat spacetimes 

To construct an asymptotically flat wormhole solution 
[23 , consider $(r) = const. Thus, from Eq. JJJ one 
obtains 



b(r) = roir/ro)- 1 ^ 



(8) 



so that b(r)/r = (r /r) (1+w) /" — * for r — > 00. Wc also 
verify that b'(r) = -(l/u)(r/r )~ (1+w)/w , so that at the 
throat the condition b'(ro) — 1/M < 1 is satisfied. 
The stress-energy tensor components are given by 
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Thus, determining the parameter oj from observational 
cosmology, assuming the existence of phantom energy, 
one may theoretically construct traversable phantom 
wormholes by considering the above-mentioned form 
function and a constant redshift function. 



into account the form function given by b(r) — vq (r/ro) a , 
with < a < 1, and using Eq. J^J and Eq. (J2J, one finds 
that the redshift function is given by 



$(r) = 



3uj + 1 

l+UJ 



In 
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(11) 



where the relationship a = —l/u> is imposed (see Ref. 
[22j for details). The stress-energy tensor components 
are provided by 



P(r) = up{r) 



3+77 



(12) 



As noted above, the spacetime is not asymptotically 
flat. Nevertheless, one may match the interior wormhole 
solution to an exterior vacuum spacetime at a finite junc- 
tion surface. 



III. STABILITY ANALYSIS 
A. Junction conditions 

We shall model specific phantom wormholes by match- 
ing an interior traversal wormhole geometry, satisfying 
the equation of state p r — cup with u < — 1, with an 
exterior Schwarzschild solution at a junction interface E, 
situated outside the event horizon, a > = 2M. 

Using the Darmois-Israel formalism 0, El , the sur- 
face stress-energy tensor, S l j, at the junction interface S 
is provided by the Lanczos equations 



1 



(13) 



where the discontinuity of the extrinsic curvatures 

across the surface E, i.e., = Kf- — K^-. The extrinsic 
curvature is defined as — n M; „ ef^e^ , where n M is the 
unit normal 4— vector to E, and e%. are the components 
of the holonomic basis vectors tangent to E. 

Taking into account the wormhole spacetime metric 
and the Schwarzschild solution, the non-trivial com- 
ponents of the extrinsic curvature are given by 
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Isotropic pressure, p r = p t — p 

It was found that considering an isotropic pressure, 
Pr = Pt — P, for $(r) to be finite one cannot construct 
asymptotically flat traversable wormholes [2^] . By taking 



and 
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The Lanczos equation, Eq. I|13[) . then provide us with 
the following expressions for the surface stresses 



V = 




6(a) 



« 2 (18) 



+ ad — 



a 2 (b-b'a) 
2(a-b) 



1 



b(a) 



,(19) 



where a and V are the surface energy density and the 
tangential surface pressure, respectively. 

We shall make use of the conservation identity, which 
is obtained from the second contracted Gauss-Kodazzi 
equation or the "ADM" constraint G^ v e^n v = Kl^ — 
K,i with the Lanczos equations, and is given by 

+ 



T p m n l 



(20) 



The momentum flux term in the right hand side corre- 
sponds to the net discontinuity in the momentum which 
impinges on the shell. 

Using S l r{i = - [cr + 2a(a + V)/a], Eq. (HJ) provides 
us with 

a' = --(a + P)+E, (21) 
a 

where H, defined for notational convenience, is given by 
b'a-b 



1 



Ana 2 



2a (1 



1 - 



(22) 



For self-completeness, we shall also include the a + V 
term, which is given by 



(J + V 
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Thus, taking into account Eq. I)23[l. and the definition of 
H, we verify that Eq. I|21|l finally takes the form 
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a^ 



which, evaluated at a static solution ao, shall play a fun- 
damental role in determining the stability regions. Note 
that Eq. I|24|) can also be deduced by taking the radial 
derivative of the surface energy density, Eq. (|T%|) . 



B. Equation of motion 

Rearranging Eq. I|18() into the form 



1 



2M 
a 



a 2 = \ 1 



6(a) 
a 



a 2 — Airaa , 



we deduce the thin shell's equation of motion, i.e., 
a 2 + V{a) = , 



with the potential given by 



2a 



M ■ 



Ha.) 



(25) 



(26) 



(27) 



where m s = Airaa 2 is the surface mass of the thin shell. 
However, for computational purposes and notational con- 
venience, we define the following factors 



F(a) = 1- 



6(a)/2 + M 



G(a) 



M-b(a)/2 



so that the potential V(a) takes the form 
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V(a) = F(a)-[ ta] - - 



(28) 
(29) 

(30) 



Linearizing around a stable solution situated at ao, we 
consider a Taylor expansion of V{a) around ao to second 
order, given by 

V(a) = V(a Q ) + V'(a )(a-a ) 

+ ±V"(a Q )(a - a ) 2 + 0[(a - a ) 3 ] . (31) 

The first and second derivatives of V{a) are given by 



V"{a) 



F 



2a 
aG\ r 



2a 

2 



- 2 



iG 



iG\[ 



(32) 



2 2a ) V 2a 



,' aG\ ( aG N , 
2|— — , 33 



respectively. Evaluated at the static solution, at a = ao, 
we verify that V(ao) = and V'(ao) = 0. From the 
condition V'(a ) = 0, one extracts the following useful 
equilibrium relationship 



m s 
2a 



m. s 



F ,_ 2 (a_oG\ /ooGy 
\m s J \ m. 



,(34) 



which will be used in determining the master equation, 
responsible for dictating the stable equilibrium configu- 
rations. 
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The solution is stable if and only if V(a) has a local 
minimum at ao and V"(a ) > is verified. The latter 
condition takes the form 



e)er<»-'' 



where VP is defined as 



(35) 



We shall now model the phantom wormhole geome- 
tries by choosing the specific form and redshift functions 
considered in Ref. [22], and consequently determine the 
stability regions dictated by the inequalities (|41l) - (|42[) . 
In the specific cases that follow, the explicit form of 
is extremely messy, so that as in |33) . we find it more 
instructive to show the stability regions graphically. 
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C. The master equation 

Using m s — Aira 2 a : and taking into account the radial 
derivative of a' , Eq. (|21|) can be rearranged to provide 
the following relationship 



(37) 



with the parameter r\ defined as r\ = V'/a', and T given 
by 



47T 

r = — (<j + v) 

a 



2tto 5' . 



(38) 



Equation (|37|l will play a fundamental role in determining 
the stability regions of the respective solutions. Note that 
the parameter ^/rj is normally interpreted as the speed of 
sound, so that one would expect that < i] < 1, based 
on the requirement that the speed of sound should not 
exceed the speed of light. However, in the presence of 
exotic matter this cannot naively de done so. Therefore, 
in this work the above range will be relaxed. We refer 
the reader to Ref. [30j for an extensive discussion on the 
respective physical interpretation of rj in the presence of 
exotic matter. 

We shall use r\ as a parametrization of the stable equi- 
librium, so that there is no need to specify a surface equa- 
tion of state. Thus, substituting Eq. (|57|) into Eq. (|33|) . 
one deduces the master equation given by 



a' m s rjo > 6 , 
where i] — 77(00) and O is defined as 



6 = — (r 2 — + — m s T . 
2tt V > Att 



(39) 



(40) 



Now, from the master equation we find that the stable 
equilibrium regions are dictated by the following inequal- 
ities 



?7o > e, 

770 < e, 



if 
if 



a' m s > , 
a' to, < . 



with the definition 
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(41) 
(42) 



(43) 



D. Stability regions 

Asymptotically flat spacetimes 

Consider the specific choices for the redshift and form 
functions given by 



$(r) = const , 
b(r) = r {r/r Q )- l l" , 



(44) 
(45) 



respectively. These are solutions to Eq. J7J), for an 
asymptotically flat spacetime. 

The factor related to the net discontinuity of the mo- 
mentum flux impinging on the shell, S, is provided by 
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The factor deduced from the equilibrium condition, T, is 
given by 
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(47) 



The radial derivative of the surface energy density, er', 
evaluated at the static solution, which will be fundamen- 
tal in determining the stability regions, takes the follow- 
ing form 




( 2n) ) (a^) 



(48) 



We shall not write down the explicit forms of the re- 
maining functions, i.e., T, \& and O, as they are extremely 
lengthy. However, the stability regions shall be shown 
graphically. 

To determine the stability regions of this solution, we 
shall separate the cases of b(ao) < 2M and b(a ) > 2M. 
From Eq. I)18|l and the definition of m s = Ana^cr, this 
corresponds to m s > and m s < 0, respectively. Here, 
we shall relax the condition that the surface energy den- 
sity be positive, as in considering traversable wormhole 
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FIG. 1: Plots for a positive surface energy density, i.e., 
b(ao) < 2M . We have defined a — ao/M, and considered 
ui = — 2 for both cases. The first plot is given by ro/M = 1, 
and the second by ro/M — 0.25. The stability regions are 
given below the solid curve. See the text for details. 
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FIG. 2: Plots for a negative surface energy density, consider- 
ing ro/M < 2. We have defined a = ao/M, and considered 
oj = —2 for both cases. The first plot is given by ro/M = 0.5, 
and the second by ro/M = 1.5. The stability regions are given 
above the first solid curve, and below the second solid curve. 
See the text for details. 



geometries, one is already dealing with exotic matter. 
Note that for a < 0, the weak energy condition is readily 
violated. 

For b(ao) < 2M, i.e., for a positive surface energy den- 
sity, and using the form function, Eq. 1)45)1 . we need to 
impose the condition rn < 2M, so that the junction ra- 
dius lies outside the event horizon, an > 2M . Thus, the 
junction radius lies in the following range 

/2M\ 

2M<a <2Ml J . (49) 

For a fixed value of u, we verify that as r — > 0, then 
an —>■ oo. The range decreases, i.e., an — * 2M, as rn — > 
2M. Note that by fixing r and decreasing u, the range 
of an is also significantly increased. 

For a fixed value of the parameter, for instance lu = —2, 
we shall consider the following cases: tq/M = 1.0, so 
that 2 < an/A/ < 4; and for r /M = 0.25, we have 
2 < ao/M < 16. The respective stability regions are de- 
picted in Fig. n From Eq. we find that a' < 0, and 
as we are considering a positive surface energy density, 
this implies m s a' < 0. Thus, the stability regions, dic- 
tated by the inequality 1)42 [I. lie beneath the solid lines 
in the plots of Fig. ^ Note that for decreasing values 
of ro/M, despite the fact that the range of an increases, 



the values of 770 & re further restricted. Thus, adopting a 
conservative point of view, using positive surface energy 
densities, we note that stable phantom wormhole geome- 
tries may be found well within the bound of < T)o < 1, 
and the stability regions increase for increasing values of 
r /M. 

For b(ao) > 2M, the surface mass of the thin shell is 
negative, m s (an) < 0. We shall separate the cases of 
r < 2M and r > 2M. 

If rn < 2M, the range of the junction radius is given 
by 

a >2M(— J . (50) 

For this specific case, a' possesses one real positive root, 
R, in the range of Eq. JSUJ), signalling the presence of 
an asymptote, g'\r = 0. We verify that a' < for 
2M(2M/r )-( 1+ ") < a < R, and a' > for a > R. 
Thus, the stability regions are given by 

— /2M\" (1-M 
r/o > 6, if 2Ml — \ <a Q <R, (51) 

rjo < 6, if a > R. (52) 



7 



10 

8~. 

6: 

4 

2 
110 
-2 
-4- 
-6 1 



co = -2 
ro/M = 1 



Stability region 




2.3 a 



FIG. 3: Plots for an isotropic pressure phantom wormhole. 
We have defined a = ao/M and considered uj = — 2 for both 
cases. For 6(ao) < 2M, the condition ro < 2M is imposed. 
The first plot is given by ro/M — 1.0, and the second by 
ro/M = 1.5. The stability regions are given below the solid 
curves. 
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FIG. 4: Plots for an isotropic pressure phantom wormhole, 
for b(ao) > 2M and r < 2M. We have defined a = a /M 
and considered w = —2 for both cases. The first plot is given 
by ro/M = 1.0, and the second by r /M = 1.5. The stability 
regions are given above the first solid curve, and below the 
second solid curve. 



Consider for uj = —2, the particular cases of ro/M = 
0.5, so that ao/M > 8, and r /M = 1.5, so that a /M > 
2.667. The asymptotes, a'\n = 0, for these cases exist 
at R/M ~ 13.9 and R/M ~ 4.24, respectively. These 
cases are represented in Fig. |5J Note that for increasing 
values of ro/M, the range of an decreases, and the values 
of 770 are less restricted. Thus, one may conclude that the 
stability regions increase, for increasing values of ro/M . 

If ro > 2M, then obviously an > ro. We verify that 
a' > 0, and consequently m s a' < 0, so that the stability 
region is given by inequality i|42|) . We verify that the 
values of 770 are always negative. However, by increasing 
ro/M, the values of 770 become less restricted, and the 
range of an decreases. 



Isotropic pressure, p r = Pt ~ p 
Consider the following functions 



*(r) = (f^) ln(r/r ) 
b(r) = r (r/r )- 1/uJ , 



(53) 
(54) 



which are solutions of a phantom wormhole possessing 
an isotropic pressure |22j. 

The factor related to the momentum flux term, 3, is 
given by 



87raQ 



1+uA 


( IQ. 


UJ > 


\ a.0 


1- 







6to 



(55) 

The r and a' are identical to the previous case of an 
asymptotically flat spacetime, and as before we shall not 
show the specific forms of the remaining functions, as 
they are extremely lengthly. 

To determine the stability regions of this solution, as in 
the previous case, we shall separate the cases of b(ao) < 
2M and b(a ) > 2M. 

For b(ao) < 2M, we have m s > 0, and the condition 
ro < 2M is imposed. Therefore, the junction radius lies 
in the same range as the previous case, i.e., Eq. I|49|) . We 
also verify that a' < in the respective range. Thus the 
stability region is given by 



r/ < 9, if 2M <a < 2M 



/2M 
\ r 



-(i-M 



(56) 
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FIG. 5: Plot for an isotropic pressure phantom wormhole, 
considering a negative surface energy density, with ro/M = 
2.5. We have defined a — ao/M and considered u = — 2. The 
stability region is given below the solid curve. See the text 
for details. 

Consider, for simplicity, oj — —2, and the cases for 
ro/M = 1 and ro/M = 1.5 are analyzed in Fig. The 
ranges are given by 2 < a /M < 4 and 2 < a /M < 
2.667, respectively. Note that as ro/M decreases, the 
range of an increases. However, the values of the param- 
eter 770 become more restricted. Thus, one may conclude 
that the stability regions increase, as ro/M increases. 

For b(ao) > 2M, then m s (an) < 0. As before, we shall 
separate the cases of rn < 2M and rn > 2M. For rn < 
2M, the range of a is given by a > 2M(2M/r )~ ( ' 1+u '\ 
as in the previous case of the asymptotically flat worm- 
hole spacetime. 

For this case a' also possesses one real positive root, 
R, in the respective range. We have a' < for 
(2M/r )~ 1/uJ < a < R, and a' > for a > R. The sta- 
bility regions are also given by the conditions (|51|) - (|52|l . 
We have considered the specific cases of ro/M = 1 so 
that the respective range is ao/M > 4; and ro/M = 1.5, 
so that ao/M > 2.667. The asymptotes, = 0, for 
these cases exist at R/M ~ 6.72 and R/M ~ 4.24, re- 
spectively. This analysis is depicted in the plots of Fig. 0] 
Note that the plots given by are inverted relatively to 
the asymptotically fiat spacetime. For decreasing values 
of r /M, note that the value of the stability parameter r/o 
becomes less restricted and the range of the junction ra- 
dius increases. Thus, one may conclude that the stability 
regions increase for decreasing values of ro/M. 

If r > 2M, then a > r . We find that a' > 0, 
which implies m s o~' < 0. Consider uj = —2, and the 
specific case of ro/M — 2.5, so that the stability region 
lies below the solid line in Fig. |3J We also verify that for 
increasing values of ro/M, the values of the parameter 
rjo become further restricted. Thus, one may conclude 
that the stability regions decrease for increasing values 
of ro/M. 



IV. SUMMARY AND DISCUSSION 

As the Universe is probably constituted of approxi- 
mately 70% of null energy condition violating phantom 
energy, this cosmic fluid may be used as a possible source 
to theoretically construct traversable wormholes. In fact, 
it was found that infinitesimal amounts of phantom en- 
ergy may support traversable wormholes (221 ) . In this 
paper, we have modelled phantom wormholes by match- 
ing an interior traversable wormhole geometry, satisfying 
the equation of state p — top with lo < — 1, to an exterior 
vacuum solution at a finite junction interface. We have 
analyzed the stability of these phantom wormholes, an 
issue of fundamental importance, to linearized perturba- 
tions around static solutions, by including the momen- 
tum flux term in the conservation identity. We have con- 
sidered two particularly interesting cases, namely, that of 
an asymptotically flat spacetime, and that of an isotropic 
pressure wormhole geometry. The latter solution is of 
particular interest, as the notion of phantom energy is 
that of a spatially homogeneous cosmic fluid, although it 
may be extended to inhomogeneous spherically symmet- 
ric spacetimes. We have separated the cases of positive 
and negative surface energy densities and found that the 
stable equilibrium regions may be significantly increased 
by strategically varying the wormhole throat. As we are 
considering exotic matter, we have relaxed the condition 
< T)o < 1, and found stability regions for phantom 
wormholes well beyond this range. There are several 
known examples of exotic 770 < behavior, namely the 
Casimir effect and the false vacuum 30] , so that one can- 
not a priori impose < 770 < 1 until a detailed micro- 
physical model of exotic matter is devised. 

As emphasized in Ref. j^, these stable phantom 
wormholes have far-reaching physical and cosmologi- 
cal implications. First, apart from being used for in- 
terstellar travel, they may be transformed into time- 
machines |l9l l2(| , consequently violating causality with 
the associated time travel paradoxes. Relative to the 
cosmological consequences, the existence of phantom en- 
ergy presents us with a natural scenario for traversal 
wormholes. It was shown by Gonzalez-Diaz 36], that 
due to the fact of the accelerated expansion of the Uni- 
verse, macroscopic wormholes could naturally be grown 
from the quantum foam. It was shown that the worm- 
hole's size increases by a factor which is proportional 
to the scale factor of the Universe, and still increases 
significantly if the cosmic expansion is driven by phan- 
tom energy |37| . However, it was also found that using 
wormholes modelled by thin shells accreting phantom en- 
ergy |38j |. the wormholes become asymptotically comov- 
ing with the cosmological background as the Big Rip is 
approached, so that the future of the universe is shown 
to be causal. 
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